A RIGOROUS REAL TIME FEYNMAN 
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Ken Loo 

Abstract. We will derive a rigorous real time propagator for the Non-relativistic 
Quantum Mechanic L 2 transition probability amplitude and for the Non-relativistic 
wave function. The propagator will be explicitly given in terms of the time evolution 
operator. The derivation will be for all self-adjoint nonvector potential Hamiltoni- 
ans. For systems with potential that carries at most a finite number of singularity 
and discontinuities, we will show that our propagator can be written in the form 
of a rigorous real time, time sliced Feynman path integral via improper Ricmann 
integrals. We will also derive the Feynman path integral in Nonstandard Analysis 
Formulation. Finally, we will compute the propagator for the harmonic oscillator us- 
ing the Nonstandard Analysis Feynman path integral formuluation; we will compute 
the propagator without using any knowledge of classical properties of the harmonic 
oscillator. 



I. Introduction. Since Feynman's invention of the path integral, much research 
have been done to make the real time Feynman path integral mathematically rig- 
orous(see [6], [9], [10], [13], [18], [19], and [20]). In physics, the real time, time sliced 
Feynman path integral is formally given by (see [3], [4], and [5]) 
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where the integral of the first equation in (1.1) is an improper Riemann integral, 
and the last line in (1.1) is the evolution operator operating on the dirac delta 
function's y variable. It is well known that mathematical rigor of (1.1) is lacking, 
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and we know that an integral over path space in real time can not be well defined 
with measure theory (see [6]). 

The problems with the objects in equation 1.1 are that we do not know if the 
improper Riemann integrals exist, we do not know if the k limit exists, and we 
do not know if the Fcynman path integral in 1.1 produces the propagator. In his 
paper(see [11] footnote 13), Feynman observed that by using wave functions, ill- 
defined oscillatory integrals can be given rigorous meaning. With this observation, 
we will reformulate equation 1.1 into rigorous mathematical objects that represent 
the propagator. 

From mathematics, we know that for some values of t, some propagators must 
be treated as distributions; the harmonic oscillator is one example(see [6] and [7]). 
Also, K t (x,xq) given in (1.1) is a function of x and x . Thus, it is natural to 
consider K t as a tempered distribution on the class of Schwartz test functions 
S (M. n x R"). The space of square integrable functions is a subset of the space of 
tempered distributions. If we consider the wave function as a distribution and take 
its inner produce with a test function, we can formally use (1.1) and get 



Equation 1.2 will be the fifth theorem in section II. 

As for the formal evolution of the delta function in (1.1), let us formally consider 
the following equation. 



where the functions K, F and G are given by equation. 2.1. If we are going to take 
the propagator as a distribution in the sense of (1.2), we might consider the limit 
in (1.3) as a distribution limit. Doing so produces the second theorem and in some 
sense the fourth theorem in section II(equations 2.2, 2.3, 2.6a, and 2.6b). 

In mathematics, there exists a rigorous formulation for a real time, time sliced 
Feynman path integral(see [7] and [8]), it reads 
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where ijj G L 2 , the integral in (1.4) is an improper Lebesgue integral with conver- 
gence taken in the L 2 topology, and the k limit in (1.4) is taken in the L 2 topology. 
Comparing (1.4) and (1.1), we see that rigorously we have a Feynman path integral 
that has all convergence taken in L 2 topology while in physics, formal improper 
Riemann integral and pointwise convergence is favored. 

What we will do is convert all convergences in L 2 topology into pointwise conver- 
gences in t. The idea is to use a wave function as a convergence factor as observed 
by Feynman. For simplicity, suppose / (x) G L 2 (R) , g (x, y) G L 2 (R x R) are such 
that they are bounded and continuous. Further, suppose that both 

(1.5) h(x)= / g(x,y)dy, 

J —a 

p (x) = lim / g (x, y) dy, 

a,b-,oo J_ a 



arc in L 2 (R) as a function of x. In (1.5), we take the integrals to be Lebesgue 
integrals and the limits arc taken independent of each other in the L 2 norm. Notice 
that for p (x) , we can interpret the integral as an improper Lebesgue integral with 
convergence in the L 2 topology. Let us denote X[-c,d\ to be the characteristic 
function on [— c, d]. Schwarz's inequality then implies 



(1.6) 



d rb 



f (x) p(x)dx- / / (a;) g (x, y) dxdy 

J — c J —a 

2 Hp -h\\ 2 + 1|/ - X[- c ,d]/ll2 \\h\\ 2 -► o. 



< 



Thus, we can write 



(1.7) / f(x)p(x)dx= lim / / / (x) g (x , y) dxdy , 

Jr a,b,c,d^oo J _ c J _ a 

where the limits are all taken independent of each other. Since / and g are bounded 
and continuous, the Lebesgue integral over [—a, b) x [— c, d] in (1.7) can be replaced 
by a Riemann integral. Since the limits are taken independent of each other, we 
can then interpret the right hand-side of (1.7) as an improper Riemann integral. 
If / and g carry singularities and discontinuities, care must be taken in the region 
of integration so that the replacement of Lebesgue integral with Riemann integrals 
can be done. 

The technique of converting L 2 limits into pointwise limits as illustrated above is 
what we will use to prove all the theorems in the next section. It is the foundation 
of this work. 

II. Results. In his paper(see [11] footnote 13), Feynman observed that by using 
wave functions, ill-dcfincd oscillatory integrals can be given rigorous meaning. With 
this observation, we will reformulate equation 1.1 into a rigorous mathematical 
object that represents the propagator. 

The goal of this paper is the following. First, we will elaborate on Feynman's 
observation and use wave functions to provide a convergence factor in the deriva- 
tion of a real time propagator that takes the form of an L 2 transition probability 
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amplitude. We will use wave functions to derive a real time, time sliced Feynman 
path integral. We will derive two Nonstandard Analysis formulations of the time 
sliced Feynman path integral. Finally we will compute the propagator of the har- 
monic oscillator using our Nonstandard Feynman path integral representation. We 
will assume that the reader is familiar with Nonstandard Analysis(see [13]-[17] and 
references within). 

The usual idea in using Nonstandard Analysis is to replace the time slice limit 
by a standard part(see [9], [13], [18], [19] and references within). We will derive a 
Nonstandard formulation that transfers the time slice limit into the nonstandard 
world and standard part is taken on infinitesimal parameters in wave functions. It 
was shown in [19] that for the harmonic oscillator, equation 1.1 can be cast into 
the language of Nonstandard Analysis where the time slice limit is replaced by a 
standard part. Further, using Nonstandard Analysis methods, one can rigorously 
compute the harmonic oscillator propagator without having prior knowledge of the 
classical path. We will follow the approach of [19] in the computations of this paper. 
In [19], we do not know if equation 1.1 is the propagator apriori; we are satisfied 
because the computation produced the correct results. In this paper, we will have 
a Feynman path integral representation that is known to produce the propagator 
and we will use it to compute the harmonic oscillator propagator. 

-ti 2 

What we will show is the following. Let H = — — A + V (x) = H + V (x) be 

essentially self-adjoint and the domain of H contains the Schwartz space of rapidly 
decreasing test functions. Denote the closure of H by H. Let t > and let 
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K(x,x ,v,l,t) = / G s (y,rj) 



exp 




The notation for K (x, Xo, f], 7, t) is that the evolution operator operates on the z 
variable while leaving xq fixed and the result is a function of xq and y\ finally, the 
y variable is integrated against G. We point out that K in (2.1) is in the form 
of an L 2 transition probability amplitude, but neither F nor G are wave functions 
since they are not normalized to 1 in the L 2 norm. Also, the form of K is similar 
to the form of the propagator given in [12] where Prugovecki provides a theory 
of stochastic Quantum Mechanics. It will be interesting to see the relationship 
between (2.1) and the stochastic propagator derived by Prugovecki. One immediate 
difference between (2.1) and Progovecki's formulation is that (2.1) stays within the 
popular representation of Quantum Mechanics where as Prugovecki uses a different 
representation(see [12]). 

The existence of K (x, x , f], 7, t) is immediate since both functions in the inte- 
grand are in L 2 . We will show that 

First Theorem. K (x, Xo, T], 7, t) is continuous as a function of(x,x ) G R 2 ™ and 
it is uniformly bounded as a function of (x,Xo) & M. 2n . 



A RIGOROUS REAL TIME FEYNMAN PATH INTEGRAL AND PROPAGATOR 5 



The kernel in the first theorem will play the role of an integral kernel in the 
following sense. 

Second Theorem Part a. Let <p, ip S L 2 (R™) Let H be essentially self-adjoint, 
then 

(2.2) f </>(£) 

lim / (f> (x) tp (xq) K (x, x a , t], 7, t) dxodx, 

where the notation for the integral in the right hand side of the equality in (2.2) 
means an improper Lebesgue integral with the convergence at infinity taken pointwise 
in t (see equation 5.10 for more details), and the limits are taken independent of each 
other and pointwise in t. 

Notice that in the above theorem, the convergence of the improper Lebesgue 
integral is pointwise in t as opposed to convergence in the L 2 topology in equation 
1.4. A pointwise convergence might provide computational advanages. 

We will not attempt to pass the limits in (2.2) inside the integral since some real 
time propagators does not exist as a function for all time and must be treated as 
distributions (see [1]). We will make a connection between K (x, x n ,n, 7, t) and the 
theory of distributions (see remark 6.2 and equation 6.4). On the other hand, we 
would be like to be able to pass the limits inside the improper Lebesgue integral 
when the kernel in the first theorem and the propagator for the evolution are well 
behaved. The next theorem provides us with that opportunity. 

Second Theorem Part b. Let (p,ip 6 L 2 {~\L 1 . Let H be essentially self-adjoint, 
then 

(2.3) f ct>(x) 

JR" 

lim / (f> (x) tp (xo) K (x, xq, n, 7, t) dx a dx, 

where the integral in the right hand side of the equality is a Lebesgue integral and 
all limits are taken independent of each other and pointwise in t. 

In the above theorem, the first theorem and the fact that the wave functions 
arc in L 1 provide us the opportunity to pass the limits inside the integral and 
produce the propagator for the evolution. We will do this for the harmonic oscillator 
Hamiltonian. Further, for the purpose of passing the limits, we will not attempt to 
generalize the wave functions to all of L 2 since integrating the propagator against 
two arbitrary L 2 wave functions in the sense of equation 1.2 is not always well 
defined; the free evolution propagator is one such example. 

The kernel in (2.1) can be explicitely represented by a time sliced Feynman path 
integral. 

Third Theorem. Let H be essentially self-adjoint, and the potential V be such 
that it has at most a finite number of discontinuities and singularities. Let 




-V(x-j) 



exp 



-itH 



(x) dx 



-itH\ 



(x) dx = 



6 A RIGOROUS REAL TIME FEYNMAN PATH INTEGRAL AND PROPAGATOR 



then 
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In (2.5), the integral is an improper Riemann integral and the k limit is taken 
pointwise in t. 

In the third theorem, there is no restriction on the type of discontinuities and 
singularities on the potential as long as the Hamiltonian is essentially self-adjoint 
and by improper Riemann integral, we mean a Riemann integral with convergence 
at infinity taken pointwise in t. Further, it is not necessary to put the restriction 
on potential in the above theorem and work with improper Riemann integrals if 
one is willing to work with improper Lebesgue integral as in part a of the second 
theorem(see remark 3.2), but for our purpose of computing the harmonic oscillator 
propagator(see section VIII), we will formulate the third theorem as above. 

The problem with formulating a real time, time sliced Feynman path integrals 
is that Fubini's theorem can not be applied due to the oscillatory nature of the 
integrand. We will see that the application of Fubini's theorem can be justified in 
the derivation of (2.5) because the functions F and G play the role of convergence 
factors as Feynman pointed out. 

The propagator is usually formulated for the wave function. If we wish to work 
with the wave function, we have the following. 



Fourth Theorem Part a. Let tp e L 2 ■ 

the following is true 
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l ), H be essentially self adjoint, then 



tp(x Q ) K (x,x ,r],^,t) dx , 



where the integral in (2.6a) is a Lebesgue integral and the limits are taken indepen- 
dent of each other in the L 2 topology. 

Part a of the fourth theorem above provides us with another way to deal with 
arbitary L 2 wave functions when equation 2.1 is considered as a distribution in R 2 ™ . 



Fourth Theorem Part b. Let ip, 

the following is true 



GL 2 



l ), H be essentially self adjoint, then 



(2.6b) 




(x) = 



ip(xo)K(x, a? ,7?,7,i) dx Q \ dx 



where the integrals in (2.6b) are iterated Lebesgue integrals and the limits are taken 
independent of each other and pointwise in t. 



We will connect the second theorem to the theory of distributions. 
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Fifth Theorem. Let S be the space of rapidly decreasing test functions. Suppose 
H = Hq + V is essentially self-adjoint, then there exists a tempered distribution 
K t (x, x ) on S (R™ x MP) such that V</> (x) ,tp(x ) e S {R n ), 



(2.7) 



<f>(x) 



-itH 



{x)dx = j K t (x, xq) 4> (x) tp (xq) dxdxo, 



where the integral in the right hand side of equation ( 2. 7) is a distribution inner 
product. 

The second theorem above is linked to the theory of tempered distributions via 
the fifth theorem. We will prove some properties of the distribution Kt(x, xo) 
given in (2.7). We have in fact gone beyond the theory of distributions in the sense 
that the second, and fourth theorem are not just true for rapidly decreasing test 
functions, they are true for a much bigger class of functions. 

Finally, the above theorems can be cast into the language of Nonstandard Analy- 
sis. The idea of using Nonstandard Analysis to formulate the Feynman path integral 
is not new (see [9], [13], [18], [19] and references within). The usual formulation is to 
replace the time slice limit with a standard part. Following that idea, the third 
theorem can be easily reformulated in the language of Nonstandard Analysis as 
follows. 

Sixth Theorem. With the notations and conditions in the third theorem, we can 
write 



(2.8) K (x , Xo,i], 7, t) = lim K k (x, x , 77, 7, t) , 

k^oo 

where Kk is given in 2.5. Let uj e *N — N, then 

(2.9) K (x, x , T), 7, t) = st {*K U (x, xo, n, 7, t)) , 

where st is the standard part. 

Notice that *K UJ consists of an infinite(namely to) copies of ^-improper Riemann 
integrals. Equation 2.9 is in fact a special case of the work done in reference [9]. 

We now come to a new formulation of the Feynman path integral in the language 
of Nonstandard Analysis. 

Seventh Theorem. Under the conditions of the second theorem(part a or part b), 
let i] and 7 be positive infinitesimals in the language of Nonstandard Analysis, then 
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where equations 2.10a and 2.10.b corresponds to theorem two part a and part 
b respectively, st means the standard part of the *-transformed improper Lebesgue 
integral(2.10a, in the sense in second theorem part a) and Lebesgue integral(2.10b, 
second theorem part b). 

The implication of the seventh theorem on Feynman path integrals is the follow- 
ing. Suppose the Hamiltonian has a finite number of singularities and discontinu- 
ities, then the third theorem holds. *-transforming the third theorem, equation 2.5 
reads: for all 77,7 g *M+, 

(2.11) *K (x, x , n, 7, t) = * lim *K k (x, x ,V, 1, t) , 

k — >oo 

where the *-limit is a limit taken in the nonstandard world, and *K k is k copies 
of *-improper Riemann integrals. In particular, we can let 77 and 7 be positive 
infinitesimals and use equation 2.11 in the seventh theorem. The time sliced Feyn- 
man path integral now has time sliced limits taken in the nonstandard world and 
standard parts taken on 77 and 7 in the sense of the seventh theorem. As mentioned 
earlier, this formulation differs from the popular usage of Nonstandard analysis 
on path integrals in that the time slice limit is not replaced by a standard part. 
This new formulation uses the fact that the functions F and G behave like delta 
functions when n and 7 are positive infinitesimals. 

Lastly, we will use equation 2.9 and theorem 2 part b to compute the harmonic 
oscillator propagator. We will compute the propagator in such a way that no 
prior knowledge of the classical path is needed. In fact, the classical part of the 
propagator naturally falls out from quantum considerations. The usual method 
to compute the harmonic oscillator with the Feynman path integral is to use the 
classical path and seperate out the classical and quantum fluctuation parts (see [4] 
and [5]). From our computational point of view, the classical mechanics part comes 
purely from quantum considerations and that goes against the grain of Feynman's 
original idea that quantum mechanics come from classical mechanics via the action 
integral in the integrand of integration over path space. 

The Hamiltonian for the harmonic oscillator is H — A H x 2 . It is well 

2m 2 

7T 

known that for < t < — , the n dimensional harmonic oscillator propagator is 

A 

given by 
(2.12) 
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= (2^) 



sin At 



g (x,x ,t) = cxp I ^ -7-^ [(x% + x 2 ) cos At - 2xx ] j 
where the g (x, xq, t) is the classical part and the h (J^j is the quantum fluctuation. 
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Given 2.12, we would expect that 

(2.13) K (x,x ,r],j,t) = h Q,*) / 9 (y,Vo,t) F (y ,x a ,j) G (y,x,r)) dydy Q . 



Notice that in 2.13, the disturbance of the functions F and G affects only the 
classical part. 

We conclude this section with a comment and a summary. The reader should 
compare the similarities and diferences between the formulations above and that 
of the notion of weak integral kernels (see [2] and references within). In particular, 
one difference is that the above kernel exists for all essentially self-adjoint Hamil- 
tonians. The main purpose of this paper is to derive a rigorous theory of real time 
propagators and real time Feynman path integrals. The propagator exists for all 
essentially self-adjoint Hamiltonians and is closely related to distributions. The 
Feynman path integral exists for potentials that carry at most a finite number of 
singularities and discontinuities, it is formulated via improper Riemann integrals, 
and it can be formulated with classical analysis and Nonstandard Analysis. Lastly, 
we use Nonstandard Analysis and compute the propagator for the harmonic oscil- 
lator without prior knowledge of the classical path. 

III. Proof of Third Theorem. We start by giving a quick proof of the third 
theorem. The third theorem is a specific case of the work work done in reference 
[9] . For the full details of the proof, we refer the reader to [9] . 

We first set some notations. Suppose V is such that it has at most a finite 
number of singularities and discontinuities. Let k G N and 1 < I < k + 1. We will 
denote the interior of the Ith box by 

(3.1) A 1 = (-a[,b[)x---x(-a l n ,b l n ), 

for positive and large a's and b's. Let K = {y\ . . .y p } be the set of discontinuous 
and singular points of V. For each y q = (y\, . . . , y^) 6 K, denote the Ith box 
centered at y q by 

(3-2) ^ = W-4r^ + ir) x - x ^--i7'^ + ^r)' 

C\ &1 Cn &n 

for positive and large c's and <fs. Let 

(3.3) CW- j(J B'j. 

For arbitrary large a's, 6's, c's, and d's, C l is a box which encloses the set K and at 
each point of K, a small box centered at that point is taken out. Associated with 
C l is a set of indices 

(3.4) = {a[,..., a l n , b[, ...,b l n , c\\ . . . , . . . , 

d 1 ' 1 d 1 ' 1 d p - 1 d p ' l \ 

We will denote by {ji} — > oo to mean 

(3.5) a\ , . . . , a n , b\, . . . , 6 n , cy , . . . , c„ : , . . . , , . . . , cP n , 

d¥,...,dy,...,d\' l ,...,<pj 
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where all indices goes to infinity independent of each other. Notice that as {ji} — > 
oo, we recover R™ a.e. from C l . We will denote by X{h} the characteristic function 
on C l . Notice that for / £ I? (MP), 



(3.6) lim X {j l }f = f a.e. 

{jl}->oo 



where the limit in (3.6) is taken in the L 2 topology. Let us write 
(3.7) L>| /(h | = C l x • • • x C'\ I < h. 

Associated with -Dj./h} is a se t of indices 

(3-8) K}=L){J«}, 

and as before, we will use the notation { J/ 1 } — > oo to mean 
(3.9) {j/} -> oo, ...,{>} -> oo, 



where the indices are taken to infinity independent of each other. Finally, we will 
denote by J Q to be Riemann or improper Riemann integration over the region O 
and J Q to be Lebesgue integration over the region O. 

Theorem 3.1. The third theorem in section II is true. 

Proof. Trotter's product formula(see [7], [8], and [10]) and Schwarz's inequality im- 
plies that 

(3.10) 




where the limit in (3.10) is taken pointwise as a function of t. To the right 
of each of the operator exp ^ — ' n (^-10), we put in the identity operator 
lim^-j j^oo X{ji } for 1 < / < A: in increasing order from right to left and the limit 

is taken in the L 2 topology. Since exp ^ f~h ^) ' 6XP ( ' mmt ^P uca ti° n 

by a characteristic function are all continuous operators, we can take all the limits 
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outside of the operators and get 
(3.11) 
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X{j k } ■ ■ 



—Sk {xk+i, ■■■,xi] 



Fx (^i)7) dxi ...dx k - 



In the last equality in (3.11), we used the integral representation of the free evolution 
operator; we emphasize that all limits in (3.11) are taken in the L 2 topology and 
are taken independent of each other. Equations (3.6), (3.10), (3.11) and Schwarz's 
inequality imply that 



(3.12) 



lim w„ jfc 

k— >oo 



exp 



-itH 



(y)dy 



r . lim X{j k+1 }Gx(xk+i,v) x 
te+i}- >oo 



— S k (x k +i,- - .,xi) 



lim w n .k / e 

lim / G s (x k +i,v) x 

{Ji }-°° Jd Ij*+i\ 



Fs (xi,j) dxi ...dx k 



dx 



fc+i 



lim w n ^k lim 



exp 



— S k (x k+ i, ■■■,x 1 ) 



Fs (^1,7) dx 1 ... dxk+i- 



In (3.12), all limits inside the integrals are taken independent of each other in the 
L 2 topology and all limits taken outside of the integral are taken pointwise in t. By 
construction, the integrand 



le 



—S k (xfc+i, . . . ,a?i) 



Fxo {xwy) 



(3.13) G s (f fe+ i,r?)exp 

is bounded and continuous on D^jk+iy Hence, the Lebesgue integral over D^jk+iy 
in the last equality of (3.12) can be replaced by a Riemann integral over D^jk+iy. 

Since the { Ji +1 } limits in (3.12) are all taken independent of each other, we can 
interpret those limits and the integral as an improper Riemann integral. □ 
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Remark 3.2. It is not necessary to use improper Riemann integrals or put the 
discontinuities and singularities restriction on the potential. We forget about the 
holes centered at each elements of K as given in (3.2) and take C l — A 1 as defined 
in (3.1). Proceeding in the same manner as the above proof, we get to (3.12). At 
this point, we do not replace the Lebesgue integral with Riemann integral since the 
integrand is not necessarily bounded and continuous over the region of integration. 
We are then left with an improper Lebesgue integral in which the convergence of 
the integral is taken pointwisc in t. 

IV. Proof of First Theorem. In this section, we prove the first theorem in 
section II. Let us denote 



(4.1) 



exp 



exp 



-itV\ 


exp ^ 


-it H Q 


kh J 


kh 


-itH 




f-itV 


kh 


J exp 


v kh 



Theorem 4.1. With our previously defined notations, we have 
(4-2) \K(x,x ,v,l,t) | < Ct,„, 7 

where Ct )J? , 7 is a constant depending only on t, n, and 7. 

Proof. Since the evolution operator has norm 1, using Schwarz's inequality on the 
kernel in 2.1 gives 



(4.3) 



\K (x,x ,v,l,t) I < 

\\G S {y,v) \\ 2 \\ F s (j/>7) || 2 = c 't.,vn 



□ 



We will now show that K (x, xo, rj, 7, t, ) is continuous as functions of (x, xq) G 



Lemma 4.2. Let f,g e L 2 (R n ), then the following is true 
(4.4) 



exp 



-it H c 
kh 



f 



(x) dx 



exp 



-itH 
kh 



(x) f (f ) dx. 



Proof. Let \a be the characteristic function of the cube centered at the origin with 
sides of length a, then 



(4.5) 



exp 



-ieHn 



(x) 



lim tu„,i / Xa (y) exp 



tme fx — y 
~2h~ 



f (y) dy, 
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where the limit in (4.5) is taken in the L 2 norm. Using Schwarz's inequality on a 
and Lebesgue's dominating convergence theorem on /?, we have 



itH 



(x) dx = 



(4.6) 

lim w nA / xp{x)g{x)\ I Xa{y) exp 
lim w nA Xa(y)f(y)\ / Xp{x) exp 

P,a^oo J Rn J Rn 



true ( x — y 

ime ( x — y 
~2h 



f(y) dyjdx = 
g{x) dx\dy, 



where the limits are taken pointwise in t. Using Schwarz's inequality on j3 and 
Lebesgue's dominating convergence theorem on a in the last expression in (4.6) 
gives (4.4). □ 

Lemma 4.3. Let f,g e L 2 (R n ), then the following is true 
(4.7) 



exp 



-itH 

h 



I 



(x) dx 



J 



exp 



-itH 
~h~ 



(x) f (x) dx. 



Proof. Intuitively, if we think of the evolution as an exponential of the Hamiltonian 
H, we can expand the exponential in powers of H and put all powers of H from the 
function / onto the function g since H is self-adjoint. This is also true for lemma 
4.2. 

Schwarz's inequality and Trotter's formula implies 



(4.6 



exp 



-itH 
H 



(x) dx — lim 



I 9{x)[T k f] 



(x) dx, 



where the limit is taken pointwise in t and T k is given in (4.1). Using Lemma 4.2, 
and X* as defined in (4.1), we obtain 



lim / g (x) [T k f] (x) dx = lim / [f k g] (x) f (x ) dx 



(4.9) 



exp 



-itH 



(x) f (f ) dx. a 



Theorem 4.4. With our previously defined notations, the expression 



(4.10) 



K (x, x , ?7, l,t)-K (y, y , 77, 7, t) 



goes to zero as (x,xq) goes to (y,yo)- 
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Proof. We first show that K (x, x , r\, 7, t) is seperately continuous in x and x , then 
jointly continuous. Schwarz's inequality implies that 

(4.11) 

1 - 1 2 

\K{x,XQ,r},^,t) - K{y,x ,r),i,t) I < 
\\G s (z,v) ~ G ff (z,r)) II x II (2,7) II = 



I r. f- \ II 2 /" TO 



cxp 



2ft 
2 exp 



77777 I y — z 



V 



dz = C 



cxp 
/ m 



2ft 



mrj / x — z 



V 



\2irhri 



cxp 



^2 



-77777 2; 



-1 2 



-mr]z 2 \ -mr, (z + y-x) \ \ 

exp I — — \dz ) = C 7it g (y, x) 



2ft rf 



2ft 



where C 7 . t is a constant independent of x a and g(y,x) is independent of Xg. 
Using Lebesgue's dominating convergence theorem on x — > y in (4.11), we get 
Cy,tg (y,x) — > 0. Using lemma 4.3 in equation 2.1, we can put the evolution oper- 
ator on Gz (y, 77). With the same reasoning as 4.11, we get 



(4.12) 



\K (x, x ,r), 7, t) -K(x, yo,V,l,t) I < 
D v .tf{yo,x ) -> 



as zo — » yo- 
Finally, 

(4.13) 



K (x, x , ??, l,t)-K (y, y , 77, 7, i) 



K (x, x , ry, l,t)-K (y, x , 77, 7, i) 



K(y, x ,r], 7,*) - K(y,y ,r), 7, i) 



< 



+ 



< 



\/ C7,*5 (y, + \/ A,,*/ (yo, xo) ->■ 



as (£,£0) -> (j7,2/b)- □ 

V. Proof of Second, and Fourth theorem. We will now prove the second and 
fourth theorem in section II. 



Proposition 5.1. Let f,g 6 L 2 , then 
(5.1) / /(*) 

JR™ 



. ^0 , 

exp — - — g 



exp 



ft 

-VH 



(z) dz = 
(z) dz = 



\ g(x)G (x, y, 77) / (y) dxdy, 
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where G{x,y,rj) is the Gaussian kernel given in equation 2.1. 
Proof. Since |/| and \g\ are in L 2 , we have 



(5.2) 



j Rn ( j Rn 15 (f } 1 |G & V) 1 lf ^ |df ) dlJ 

/ | fl (f)' L -^ 0, 



cxp 



(x) dx < oo. 



Equation 5.1 then follows from Fubini's theorem. □ 

Theorem 5.2. Part b of the second theorem in section II is true. 

Proof. First recall that in part b of the second theorem, the wave functions 4>,ip & 
L 2 n L 1 . Using lemma 4.3, the kernel in (2.1) can be written as 



(5.3) K (x,x ,r],^,t) = 
We have that 



exp 



exp 



-itH 

h 



Gx (z, rj) 



{So). 



(5.4) 



exp 



lim 



h 

exp 



itH\ 



(x) 



-VH 



exp 



-itH 



exp 



-j H 



where the limits in 5.4 are taken in L 2 topology. Using Schwarz's inequality, the 
following holds 



(5.5) 



4>{x) 

lim / 



exp 

<j>(x) 



-itH 

h 

exp 



(x) dx = 



-vHo 



cxp 



-itH 



cxp 



-jH 



(x) dx, 



where the limits in 5.5 are taken pointwise in t. Lemma 4.3 and proposition 5.1 
implies that 



(5.6) 



cxp 



-VHo 



h 

G (x, y, n) 



exp 



cxp 



-itH\ 
-itH 



J 





'-itH\ 


exp ^ 


. » ) 






exp ^ 





n 

Gx (z, rj) 



exp 



exp 
exp 

(y) 

itH 



-jH 

h 

-jH 



exp 



1> 
iHo 



(x) = 
(y)dy 



h 



1> 



(y) dy 



Gs {z, rj) 



(x ) ip (x ) dx = 



J 



K (x, x , rj, 7, t) ip (x ) dx . 
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Thus, 
(5.7) 



/ (x) exp ( — ^— J ip (x) dx = 

lim / 4>{x) / K (x,xo,r),7,t)ip(x )dx )dx 
v,-r->oj Rn \ Jr» / 



lim 



^ (x) if (x, f , ??, 7) V> (^o) dxodx, 



where the last equality in 5.7 is obtained from the first theorem in section II(theorem 
4.1 and 4.2) and the fact that the wave functions <fi and ip are in L 1 . □ 

Theorem 5.3. Part a of the second theorem in section II is true. 

Proof. Let C 1 = A 1 , C 2 — A 2 be as described in equation 3.1 and remark 3.2. Let 
X{j 1 }jX{j 2 } be the characteristic function on the region C 1 , C 2 respectively({ji}, 
and {j'2} are as described in equation 3.4 for C 1 and C 2 ). Let <fi, ip E L 2 , then 



(5.8) 



exp 



lim 

{j"i}^oc 



-r)H { 



h 

exp 



exp 

h 



-itH 
exp 



exp 
-itH 



h 



-jHo 

h 

exp 



7#o \ 1 



lim - 2 



X{ji}< 



{j 2 }->oo 

where all limits in 5.8 are taken in the L 2 topology. Using 5.4, 5.5 and 5.8, we have 
(5.9) 

-itH y 



exp 

4>(x) 



h 

exp 



-VH 



f 0(5) 
lim / 

lim / lim X{ Jl} ^( f ) x 



(x) dx = 



exp 



-itH 



lim 

{j 2 }->oo 



exp 



ft 



exp 



exp 



exp 



h 



(x) 



(x) dx. 



Finally using Schwarz's inequality and theorem 5.2 on the last expression in 5.9 
gives 
(5.10) 



/ 0(5) 

JR™ 



exp 



-itH 



(x) dx = 



lim lim 

')>7^0{j 1 },{j 2 }^oo 



exp 



-VH 



exp 



exp 



h 



(x) dx 



lim lim 

^7^0{j 1 },{j 2 }^oo t / c i xC 



> (x) if (x, xq, 7; i> (xo) dx dx 



lim 



' (x) K (x, x , »j, 7, i) V" (zo) dxodx, 
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where all limits outside of integrals are taken pointwise in t and J by definition 
is an improper Lebesgue integral with convergence at infinity taken pointwise in 
t. □ 

Remark 5.1. If we restrict <p, tp G L 2 to have a finite number of singularities and 
discontinuities as in the proof of the third theorem (theorem 3.1), we can obtain 
an improper Riemann integral as opposed to an improper Lebesgue integral for 
theorem 5.3. 

Theorem 5.4. Part a of the fourth theorem, in section II is true. 

Proof. Equation 5.6 is true for all tp G L 2 . Taken lim, 7j7 ^ 00 on both sides of 5.6 in 
the L 2 topology gives part a of theorem four in section I equation 2.6a. □ 

Theorem 5.6. Part b of the fourth theorem in section II is true. 

Proof. Let <p,tp G L 2 , then Schwarz's inequality and part a of the fourth theo- 
rem(theorem 5.4) implies 



(5.11) 



/ 

J 



<p (x) exp 



-itH 



(x) = 



4>(x)\ lim / tp (xq) K (x, xo, rj, 7, i) dx a \dx = 
^.7— oy Rn ' J 

lim / 4>(x)[ ^(x )K(x,xo,r),i,t) dx ]dx, 

where limits inside the integral are taken in L 2 and limits taken outside the integrals 
are taken pointwise in t. □ 

VI. Proof Fifth Theorem. In this section, we prove the fifth theorem in section 
II and derive some properties of the tempered distribution in the fifth theorem. 
Since we will be working with tempered distributions, we will let <j) an d ip be in 
the class of rapidly decreasing test functions which we will denote by S (R n ). If 
and tp are elements of L 2 (M n ), we can choose a sequence of test functions {(pi} and 
{tpj} such that (j>i — > <p and tpj tp in L 2 . Applying Schwarz's inequality and using 
the fact that the evolution operator has operator norm equal to 1, we can write 
(6.1) 



/ 

JR' 



exp 



-itH 



tp 



(x) dx 



lim 



(pi (x) 



exp 



-itH 



(x) dx. 



Thus, by taking limits in the sense of (6.1), we can always recover all of the L 2 
wave functions in the theory. 

Theorem 6.1. The fifth theorem in section I is true. 

Proof. Suppose -j> fc (x)} C S (W l ) with <p k (x) (x) in S (R n ), then 



(6.2) 



[<P (x) - (f>k (x) 
- <Pk H2 x II exp 



exp 



-itH 



-itH 



tp 



(x) dx 



< 



18 A RIGOROUS REAL TIME FEYNMAN PATH INTEGRAL AND PROPAGATOR 



Suppose {x)}cS (R™) with ip k (x) -> V> (x) in S (R n ), then 
(6.3) 



0(f) 



cxp I — ^— ) (^ - V'fc) 



|2 x || exp 



h 

-itH 



(x) dx 



< 



(tp - tp k ) ||2 



| 2 x 1 1 (V — V'/O II2 -> 0. 



Hence the theorem follows from Schwartz's Kernel Theorem. □ 
Remark 6.2. Notice that the second and the fifth theorem imply that 

(6.4) 

lim / <f> (x) tp (x*o) K (x , xo, rj, 7, t) dxodx = / K t (x, xo) <f> (x) ip (xo) dxdxo, 

At t = 0, the evolution operator becomes the identity operator. In distributions 
language, We have the following 

Theorem 6.3. At t = Q, K t (x, Xo) satisfies K (x, Xo) = 5 (x — Xo) 

'-itH s 



Proof. Since 



exp 



h 



1> 



(x ) — ip (x) when t — 0, we have 



(6.5) / (x, xo) <p {%) "0 {x*o) dxodx = / <j> (x) ip (x) dx = 

J JR™ 

J J (j) (x) 5 (x — xo) ip (xo) dx dx . 

We extend (6.5) to all of S (R™ x R n ). Let 7]{x,x ) G S (R n x R"). Choose a 
sequence of functions r\ k (x, x ) in S (R n x R") such that rj k — ► rj in S (R™ x R") 
and T] k = J2lto u hk (x)vi,k (xo) where u itk and v iyk G D (R n ), the C°° compactly 
supported test functions, then 

(6.6) 

/K (x, xo) r] (x, xo) dxdx — lim / K (x, xo) r) k (x, x*o) dxdxo = 
k^oo J 

lim / Ui k (x) Vi k (x) dx — / r\ (x, x) dx = / 5 (x — x*o) rj (x, xq) dxdxo- □ 

jr" j^i ' Jr" J 



i=0 



It is well known that the free propagator satisfies K\ rcc (x, Xo) = K^ cc (x , x ) . 
We will show a similar property for the tempered distribution in theorem 6.1. 
Intuitively, it is reasonable to believe that the from lemma 4.5 we can conclude the 
following. 

Theorem 6.4. K t (x,Xo) = K t (xo,x) where equality is in the sense of distribu- 
tions. 
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Proof. V0 (x) , ip (x ) £ S (IT), we have that 
(6.7) 



K t (x, xq) <f> (x) ip (xq) dxdxa — \ 4>{x) 
itH 



exp 



-0 (x ) 



exp 



(x ) dx = 



L * (a?) 



-itH 
h 



exp 



(x) dx 



-itH 



(x) dx 



(x, x ) ^ (x) (x ) dxdxo = J K t (x , x) V> (x*o) (x) 



dxdxo ■ 



We extend (6.7) to all of S (R n x R"). Let r/(x,x ) £ 5 (K n x M"). Choose a 
sequence of functions 77^ (x, Xo) in S (M 2 ™) such that r\k — > 77 in 5 (M 2 ™) and 77^ = 

Si=o Ui > fc (%) v i,k (x*o) where and u^fc G £) (R™), the C°° compactly supported 
test functions. We then have 



(6.8) 



K t (x, x ) 77 (x, x ) rfxdx = lim / K t (x, x ) 77*1 (x, x ) dxdx = 



— >oo 



lim y K t (xo 7 x)i]k(x,x )dxdx ~ J K t (x , x) 77 (x, x ) dxdx n . □ 

VII Proof of Sixth and Seventh Theorem. 

Theorem 7.1. TTie sixt/i theorem in section II is true. 

Proof. The nonstandard equivalent of equation 1.8 is: for all u> £ *N — N, 



(7.1) 



K (x, x , T), 7, t) = st (*K U (x, x , 77, 7, *)) • 1=1 



Theorem 7.2. T/ie seventh theorem in section II is true. 
Proof. Let 



(7.2) i? (r?, 7) = / (x) ^ (x ) K (x, x , ??, 7, t) dx dx, 

H (v,l)= / VK^o) A (x,x ,r], 7, t) rfx (ix, 



c 



(x) dx. 



The second theorem of section II implies that for all e £ K + , there exists a <5 6 
such that 



(7.3) 



\H — C\ < e when n,j<5. 



We now *-transform 7.3 and conclude that any positive infinitesimal 77 and 7 is less 
than S. Since we can do this for any standard e, the seventh theorem hold for H. 
A similar argument shows that the theorem is also true for H □ 
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VIII The Harmonic Oscillator. We now compute the harmonic oscillator prop- 
71" 

agator for < t < — using the formulas above. Some of the techniques that we will 
A 

use was previously worked out in [19], for full details, we will occassionally refer 
the reader to [19]. For the harmonic oscillator, equation 1.5 rcads(with a shift in 
the indices) 



(8.1) K{q,q ,r),7,t) = 



lim w n , k+1 / Fg- (x 0,7) Gg- rj) x 

fe^OO JrR(h + 2)n 



exp 



. fe+i 

It \ 

3 = 1 



rn(x, i 



in 
~2 



^A 2 {Stf 



dx . . . dxk+i- 



Let us write Xj — (a;j, . . . a;™), and 



(8.2) 



exp 



it 



fe+i 

E 




The popular method to compute the time sliced harmonic oscillator path integral 
is to use (8.2) to decouple the integrals in 1.0 and reduce the problem to produces 
of one dimensional harmonic oscillators. Due to the extra xq, Xk+i integrals in 8.1, 
it is not immediately clear that we can use 8.2 to decouple the improper Riemann 
integrals. 

For the moment, let us consider just one of the entries in the product of 8.2. To 
shorten notation, let us write 



(8.3) 



. fc+i 

It X 

im 
2ht 

k 



(x^) 2xqX± + (x^ +1 ) 2x%x'j* +1 - 



k fe+i 

^2(,«) 2 ! -E 2 4^-i-^ 2 E(4) 
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/o ... 
1 



2 \2 



e z X 



( 1 

-1 




V o 

. . 



-1 




... 

Vo ... 

/o 

2-10 ... 
-12-10 
0-12-1 





. 
. 

Vo . 



... 0\ 



-1 

-1 1/ 

. ... 0\ 
... 



1 
1/ 









-12-10 
0-120 
0/ 



where Tk is the (fc + 2) by (fc + 2) symmetric matrix, 



(8.4) 



/ 1 

-1 


V o 



-1 



Sk 



0\ 


-1 

-1 l-e 2 A 2 / 



with Sk being the k by k symmetric matrix Sk = Ak — e X Bk, where 



(8.5) 



/ 2 -1 ... 
-1 2 -1 
0-12-1 






V o 



... o x 

... 
... 







0-12-1 
... 0-12/ 
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/l 
1 



Bi. — 









Vo 



and x a is the column vector 



1 
1/ 



(8.6) 



Let w a : [0,t] -> R be such that |w a (s) | < oo and w a (0) = x%,w a (t) = x% +1 . In 
the literature, the path w a is usually taken to be the path of the classical harmonic 
oscillator. Here, we allow w a to be any finite path that starts at Xq and ends at 
x k+i- We do not assume prior knowledge of classical mechanics. We make the 



substitution xf = 



,J " (l+l) + yf = w f + y? (notice that y% = 



since 



w a (0) = Xq and w a (t) = x% +1 ). Using the fact that T k is symmetric, we have 



(8.7) 



(xTf T k x° = (r + w a f n (r + w a ) = 

(w a f r k w a + (rf nr + («r)' nr + arY nw a = 
(w a ) l T k w a + (rfnr + ewf ^ + (ki'w)' 
(w a f r k w a + (rf nr + {T k w a f r + m* nr = 



where 



(8.8) 



\ / 



w 



By using yft = Q = y k+1 and writing T k as 

... 0\ 



(8.9) 



/ 


s fe 
V ... 






0/ 



u k+l 



( 1 

-1 



V 



^0 \ 



L k+1 



-1 




0\ 


-1 

-1 l-e 2 A 2 / 
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we obtain 
(8.10) 

where 
(8.11) 



(w a f T k w a + (y a )* S k y a + 2 y a 









i < \ 




y a = 


\yt) 




\ w% J 





= S k w a — w 



Lemma 8.1. Let f e I and < t < — . For any uj e *N - N, *S U is positive 

A 

definite in the * -transformed sense. Here, *S UJ is the *-transform of the matrix S k 
defined after equation 8.4- 

Proof. Sec [19]. 

We now go back to equations 8.1 and 8.2 in Nonstandard Analysis form. With 
an abuse of notation, in Nonstandard Analysis 8.1 reads 



(8.12) 



K{q,qQ,r],'),t) 



st\w n ^ +1 F ?0 (x ,^)G ? (x OJ+1 ,r]) 

( ■ w+1 

exp !E 



m I Xj — Xj-i 
~2 



171 \2 /- \2 



dx ■ ■ ■ dx^+i 



rtl w„, w+ i / Ffo (x , 7)Gg-(f u;+ i, rf) x 



m / Xa — Xj-i \ m 



dx . . . dx u+ i > = 



f[ exp | (^-) ((£*)< T^) \dx ...dx, 

a— 1 k \ / ^ 



We perform a ^-transform of the change of variable described in equation 8.7 on 
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and obtain 

(8.13) 

st I w niU+ i / Ff (x , 7)G f (x a ,+i, rf) x 



a=l 



a=l 



n exp i ( We i \ {r)t Su)r + 2 (pQ)t r l r d ^ d ^+^ ■ ■ ■ d $< 



Since 5 W is positive definite, it is intertable. Since S u is symmetric, the following 
is true 

(8.14) (y a f S u y a + 2 (/?*)* = (y« + Sj 1 /^)' 5 W (y Q + S^fT) - (fTf S^fT. 

Using 8.14 in 8.13 and performing the transformation z" = y" + (S^ 1 p") , we 
obtain 



(8.15) 



K {q,qo,V,l,t) 
St 



n 

il eXP 



a=l 



( im 
\2he 

(im \ 
We) 



(i* 1 )* S^z"- > dxodx^+idzi . . . dz^ 



Notice that before the improper limits are taken on the integrals, the limits of 
integration on the variables zf are dependent on x an d x u+ \ due to the fact that 
(P is dependent on them. Thus, we still can not decouple the improper Riemann 
integrals. 

Let us take a look at the limits of integrations more closely. Before the improper 
limits on the integrals are taken in equation 8.15, we have 

(8.16) w n , u+1 / F iro (xo,'y)Gj i (x u+1 ,T])x 
Jo 



a=l 



II- exp I ( We ' S "^* I d ^ ' ' ' d * w \ dS ° d ^+ 1 ' 



where both O and O are *-compact and the boundry of O depends on x , x u+ i 
and a set of indices {J} such that as { J} — > oo, O — > *K wn in the ^-transformed 
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sense. The reason for O's dependent on xq , x^+i is due to the fact that p?* is 
dependent on them (equation 8.11) and we performed the change of variables from 
equation 8.13 to equation 8.15. Further, the boundary of O is also indexed by a 
set similar to that of {J}. What we would like to do is pass the {J} limits inside 
the O integral and decouple the improper Riemann integrals in equation 8.15 into 
improper Riemann integrals in dz\ . . . dz u then dx*odx w+ \. 

It is well known from the one dimensional harmonic oscillator that 



Let us fix an O. Since O is compact (see the construction in equation 3.12). We 
can *-transform and conclude from 8.16 and 8.17 that for any (3 e *M+, there exists 
a fixed M e *M + that depends only on O such that 



whenever all entries of {</} are bigger than M. In other words, because O is 
compact, for all (xo,%w+i) 6 O, equation 8.18 is true whenever all entries of {J} 
are bigger than a fixed M; further, this M depends on O. 

Equation 8.18 allows us to use *-Lebesgue dominating convergence theorem and 
pass the {J} limits inside the O integral and decouples the improper Riemann 
integrals. Thus, we have proved the following 



(8.17) 




(8.18) 
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Theorem 8.2. For the harmonic oscillator, 
(8.19) 

K (q,qo,r/,j,t) = 

st\w„ tU!+1 F fo (x , i)G<i{x u +i,r}) 



| ( im 
2he 



(w a ) z T u w a -(fP f S-'ff 



'n,w+lS / 



st< w 



F qo (^0,7) ^(^+1,7?) X 



Q = l 



a=l 

m 



n e ^ 1 ( H ) 5 ^ r dfi • • • 



2-Kihe 



det S u 



uj JrR 2 



F qo (^0,7) G$(X U+1 ,T]) X 



2he 



(w a ) z T u w a -((?*)* S^jF 



dx x { 



w+l 



Proof. See above. □ 

It remains now to compute the last equality in 8.19. 
Proposition 8.3. With the previously defined notations, 



(8.20) 



st ■ 



/ m \ 2 
V 271-^6 7 



1 



/ m \ f / A 
\2irihJ \ sin At 



Proof. See [19]. □ 

Proposition 8.4. Let Xo,x n = y G R n &e /iied. VFii/i i/ie previously defined 
notations, 

(8.21) 



k a— 1 ^ x 7 

^ { T sTmAi + ^ C ° S A< ~ 2 ^°] } ' 

Remark 8.5. Notice that proposition 8.4 does not take place in the nonstandard 
world. 

Proof. This is just the classical version of the nonstandard results obtained from 
[19]. See [19] for more details. 
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Proposition 8.5. With the previously defined notations, 
(8.22) 

[ JrR 2 ™ 



II ex P { \2f-) {^f^uT - {(Pf S-^j dx x^ +1 j = 
(f ,7) Gg(y,v) x 

JR 2 ™ 

exp | ^ . ^ [(xq + y 2 ) cos At - 2yx ] \ dx a dy. 
^ sin At J 

Proof. Using proposition 8.4 and Lebesgue's dominating convergence theorem, we 
obtain 

(8.23) 

st \ / F q (zo, 7)G f (x a ,+i, 7?) X 

[ JrR 2 " 

ftexp {(£) [(^'T^-^^V- 

a=l ^ ' 

lim / _Fg- 

JR 2 " 

a=l ^ ' 

JR 2 ™ 

lim { I] cxp { f ^ 7W« - (/?*)* S^/F"] j ) 



dxody = 



2& 



/* r iTTi -\ i 

/ ^Vo (^o , 7) Gq (y, V) exp <^ — ^— [ (f o + V 2 ) cos At - 2yx ] \ dx dy. □ 
j R 2n ^ ft sin At J 



Theorem 8.6. for the harmonic oscillator, 
(8.24) 



2irih/ V sin At 



/ m 
\2riJ 

J 

JR 2 " 



A \ 2 



^Vo (^o, 7) Cg- (y, 77) exp <{ -Y^j-j: [(^o + J/ 2 ) cos A * - 2 y f o] }> dfbdy. 



Proof. Equation 8.24 follows from theorem 8.2, proposition 8.3, and proposition 
8.5. □ 
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Theorem 8.7. Let <p,tjj E L 1 n L 2 , then for the harmonic oscillator Hamiltonian 
and for < t < n/X, 



Proof. Notice that \K (q,q ,r],j,t) \ < C n ,\,t- Substituting equation 8.24 in equa- 
tion 1.3 and using Lebesgue dominating theorem on the 77,7 limits give 8.25. □ 
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